ASP Conference Series, Vol. **VOLUME***, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 

The physics of extrasolar gaseous planets : from theory to 
observable signatures 

G. Chabrier, F. Allard, I. Baraffe 

C.R.A.L., Ecole Normale Superieure de Lyon, Lyon, France 
T. Barman 

Wichita State University 
P.H. Hauschildt 

Hamburger Sternwarte, Hamburg 

Abstract. We review our present understanding of the physical prop- 
erties of substellar objects, brown dwarfs and irradiated or non-irradiated 
gaseous exoplanets. This includes a description of their internal proper- 
ties, mechanical structure and heat content, their atmospheric properties, 
thermal profile and emergent spectrum, and their evolution, in particular 
as irradiated companions of a close parent star. The general theory can 
be used to make predictions in term of detectability for the future obser- 
vational projects. Special attention is devoted to the evolution of the two 
presently detected transit planets, HD209458B and OGLE-TR-56B. For 
this latter, we present a consistent evolution for its recently revised mass 
and show that we reproduce the observed radius within its error bars. 
We briefly discuss differences between brown dwarfs and gaseous planets, 
both in terms of mass function and formation process. We outline sev- 
eral arguments to show that the minimum mass for deuterium burning, 
recently adopted officially as the limit to distinguish the two types of ob- 
jects, is unlikely to play any specific role in star formation, so that such 
a limit is of purely semantic nature and is not supported by a physical 
justification. 



1. Introduction 

Substellar objects (SSO) are defined as objects not massive enough to sustain 
nuclear fusion in their core, i.e. with mass below the hydrogen-burning minimum 
mass {ruHBMM ~ 0.072 M©, Chabrier &: Baraffe 1997; see Chabrier & Baraffe 
2000, Burrows et al. 2001 for recent reviews). These once elusive objects have 
been revealed to be about as numerous as the Galactic stellar population. Several 
hundreds of brown dwarfs (BD) have now been discovered in the field or in young 
clusters and more than a hundred gaseous planets (CP) have been detected 
outside our solar system. Many ground-based and space-based searches are 
planed for the next coming years to seek more of these objects, with as a holy 
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grail the direct observation of an extrasolar planet. These projects require a 
sound theory to understand the physics of SSOs and to provide reliable guidance 
to observation. In this paper, we examine the present status of this theory and 
confront it to available observations. 

2. Equation of state. Internal structure 

Internal conditions for SSOs, from the HBMM to jovian conditions, range from 
Tc ^ 10^ K to 10"^ K and from pc ~ 10^ g cm~^ to 10 g cm~^ for the central 
temperature and density, respectively. Under these conditions, the average po- 
tential energy {Ze)'^ /ai of two ions separated by a distance o, — 'W/'inNi can be 
up to 50 times larger than the average kinetic energy kT. This is illustrated by 
a so-called plasma coupling parameter T = {Zef' j aikT > 1, meaning that the 
ionic plasma strongly departs from the perfect gas conditions encountered inside 
massive stars. As a comparison, the interior of Sun-like stars is characterized 
by r ~ 0.1. Under the aforementioned interior conditions, the Fermi energy of 
the electrons is comparable to the thermal energy of the plasma, kTp kT, so 
that the electrons form a finitc-tcmpcraturc, partially degenerate plasma. This 
is illustrated by the so-called degeneracy parameter 9 = kT/kTp ~ 0.05-1. This 
differs noticeably from for instance the case of white dwarfs, where the electrons 
arc fully degenerate, with 6 0.01. Last but not least, the average binding 
energy of an electron around a nucleus, Ze^/ai, is of the order of the electron 
Fermi energy ep, so that both H and He atoms experience pressure- ionization, 
a long standing problem in dense matter physics. The correct description of the 
thermodynamic properties of the strongly correlated, finite-temperature, quan- 
tum electron-ion plasma which composes the interior of SSOs thus requires the 
solution of the so-called N-body problem, one fundamental problem of statistical 
physics. Until recently, this challenging problem could not be probed by experi- 
mental conditions on Earth. Things changed within the past few years with the 
realization of laser-driven shock wave experiments at Livermore (Collins et al. 
1998) and Z-pinch experiments at Sandia (Knudson et al. 2001, 2003) which now 
reach the megabar range at temperatures characteristic of SSOs, providing strin- 
gent constraints on the internal equation of state (EOS) of these objects. The 
EOS determines not only the mechanical properties of the object (mass-radius 
relationship and contraction work) but also its internal heat capacity and adia- 
batic gradient. The EOS calculated a decade ago by Saumon &; Chabrier (1992, 
1995) has been shown to adequately reproduce the experimental results up to 
about 0.5 Mbar (Saumon et al. 2000). In the very high density limit, -P>10 
Mbar, this EOS adequately reproduces Monte Carlo simulations of dense cou- 
pled plasmas or first principle expansion schemes (see Saumon & Chabrier 1992). 
In the ~ 0.5-5 Mbar range, the laser-driven and Z-pinch experiments yied dif- 
ferent results, so that more experimental data are needed before the behaviour 
of hydrogen in this pressure range can be robustly determined. 

3. Atmospheric structure. Emergent spectrum 

Several recent reviews include a description of the spectral evolution of low-mass 
objects, from the Sun to jovian planet conditions (Allard et al. 1997, Chabrier & 
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Baraffe 2000, Burrows et al. 2001). The atmospheric physical properties of SSOs 
can be captured from the equihbrium condition of an hydrostatic atmosphere : 

dP 

dr = (pkr) dz = (1) 

With Rr ^ 0.01-1 cm^g-i and R ~ 0.1 R©, i.e. g = Gm/R^ « 0.1-10 g©, 
this yields Pph ^ g/R ^ 0.1-10 bar, ~ 10~^-10~^ gcm~^ at the photosphere, 
quite dense conditions for stellar atmosphere standards. As a function of Teg, 
the evolution of the spectral energy distribution can be outlined as follows: 

- Tcfj <^ 4000 K: H2 molecules form. Under the aforementioned density con- 
ditions, collisions either between H2 molecules or between H2 and He atoms yield 
induced dipolar transitions and thus a frequency-dependent, so-called collision- 
induced absorption (CIA) opacity which scales within first order as the square 
of the H2 abundance, «;c'i'A(H2-H2)oc n'\j^HH2- Besides H2, other molecules form, 
such as hydrides, FeH, MgH, CaH, metal oxydes, TiO, VO, the dominant sources 
of absorption in the optical, H2O, CO, which determine the shape of the spec- 
trum in the infrared. 

- Teff ^ 2500 K: various refractory compounds (grains) form, such as 
Mg2Si04, CaTiOs, MgAl204,..., sequestering the metals. These grains yield 
grey-like absorption in the optical and thus, by energy conservation, lead to 
a backwarming of the inner atmosphere. This destroys the aforementioned 
molecules absorbing in the IR, yielding red colors. This lack of metal oxyde 
(TiO, VO) absorption in the spectrum defines the so-called L-dwarf spectral 
class. On the other hand, neutral alkali lines (K, Li, Cs, Rb, Cs) in the optical 
and near-IR become increasingly important. 

- Teff ^ 1800 K. In this regime, the dominant equilibrium form of carbon 
becomes CH4, at the expense of CO. The absorption of methane in the H and K 
bands yields a redistribution of the flux at shorter wavelengths, yielding bluer 
colors. This presence of methane in the spectrum defines the T-dwarf spectral 
type, which encompasses both cool BDs and GPs. The benchmark of this type 
of objects, with jovian-like atmospheres, is the first unambiguously identified 
BD, G1229B (Oppenheimer et al. 1995), of which spectral energy distribution 
was successfully predicted by the theory (Allard et al. 1996, Marley et al. 1996). 

- Teff ^ 500 K. We enter the conditions characteristic of the jovian planets 
of our solar system, reaching the condensation limits of H2O and then NH3. 

Figure 1 displays the spectrum of a low-mass M-dwarf with Teg = 2500 K, 
a L-type BD with Teg = 1800 K and a T-dwarf BD with T^s = 1000 K, typical 
of non-irradiated CP conditions, from Allard et al. (2001) atmosphere models. 
This series of spectra illustrates the aforedescribed temperature evolution of the 
dominant sources of absorption. The role of grains in the L-dwarf spectrum is 
obvious from the grey-like behaviour at short wavelengths, due dominantly to 
Mie scattering, and the reduced bands of H2O and CO in the IR, compared with 
the hotter spectrum. The constant replenishment of these grains in the atmo- 
sphere is due to convcctive motions which reach the regions near the photosphere 
(Ludwig, Allard &: Hauschildt 2002). For cooler temperatures and luminosities, 
characteristic of T-dwarfs and planets, most grains form below the photosphere, 
convection retreats to deeper layers and the largest grains sink under the action 
of gravity to deeper, hotter regions where they melt and sublimate. 
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Figure 1. Spectral energy distribution of a typical M-dwarf (Tefj = 
2500 K), L-dwarf {T^s = 1800 K) and T-dwarf {T^s = 1000 K), from 
top to bottom, from Allard et al. (2001). Surface gravities and radii 
are from Chabrier et al. (2000a). 

4. Evolution 

Given the strong frequency dependence of the emergent spectrum of SSOs, and 
possible convective instability in the atmosphere, correct evolutionary calcula- 
tions require a consistent (non-grey) treatment between the internal structure 
profile and the atmospheric profile, as described in Chabrier & Baraffe (1997). 
A grid of atmosphere models with various (log(7, Tgff) conditions is computed, 
and matched to the interior profile at deep enough optical depth (r ^ 100) to 
reach the internal adiabat. For a given mass and age {m,t), there is a unique 
solution matching the two profiles, yielding consistency between the observable 
signature, spectrum and magnitude, and the fundamental properties of the ob- 
ject, m, R, Tcfj, and t. The evolution obeys the usual first and second principles 
of thermodynamics: 



where S denotes the specific entropy of the BD or the planet and R its radius. 

The evolution of the mass-radius relationship in the stellar and substellar 
regime is shown in Figure 3 of Chabrier & Baraffe (2000). The relatively con- 
stant behaviour of the radius as a function of mass in the substellar regime stems 
from the compensating contributions of the quantum electron gas, which yields 
an increasing radius with decreasing mass, a well-known property of degener- 
ate matter, and of the classical, strongly correlated ionic plasma. This yields a 
maximum radius near 4 Mjup (Hubbard 1994, Chabrier & Baraffe 2000). This 



L = AnR^aT^ff = - 



f du P dp 
J ^dt dt' 



f r^dS , 
J dt 



(2) 
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m-R relation has been confirmed recently in the vcry-low-mass star range by the 
remarkable agreement with the radii measured by interferometry (Lane, Boden 
Sz Kulkarni 2001, Segransan et al. 2003). On the other hand, recent observa- 
tions of the radii of the eclipsing binary YY-Gem, with masses mi=m2=O.6M0 
(Torres Sz Ribas 2002) show a ~ 10% discrepancy between the theoretical and 
observed radii. These stars, however, show intense signatures of activity, a com- 
mon feature for very short period systems such as eclipsing binaries. The surface, 
indeed, is observed to be covered by magnetic spots, reducing the total irradiat- 
ing area. Calculations in the case of cataclysmic variables (Rcnvoize et al. 2002) 
have shown that a reduced radiating surface leads to a reduced contraction at 
early stages of evolution, and thus a larger radius. Similar effects may apply to 
active, short-period systems such as YY Gem. 

Evolutionary sequences L{t), Tes{t) and color-magnitude diagrams for SSOs 
down to 0.5 X 1O~^M0 = 0.5Mj„p are presented in Baraffe et al. (2003). In 
a CMD, the Teff -sequence of absorbers described in §3 leads to red ward and 
blueward zig-zag behaviours along an isochrone. 

5. Effect of irradiation 

The majority of the presently observed exoplanets are in very close orbit to their 
parent star. In that case, the incoming stellar radiation modifies substantially 
not only the emergent spectrum but also the evolution of the planet, as shown 
initially by Guillot et al. (1996). These modifications are examined below. 

5.1. Atmosphere 

Under strong irradiation conditions, the solution of the transfer equation, which 
determines the atmospheric thermal structure of the planet, must include the 
incoming incident stellar flux in the source function. Such calculations were first 
conducted by Seager & Sasselov (1998), using a limited set of opacities, and have 
been improved by Barman et al. (2001) and Sudarsky et al. (2003). Figure 5 of 
Baraffe et al. (2003) compares the thermal structure of two planets with intrinsic 
temperatures T^g = 1000 K and T^fi = 100 K, respectively, orbiting at 0.05 AU 
of a G2 star, when irradiation is or is not taken into account, as calculated in 
Barman et al. (2001). The irradiated thermal profile is shown to be strongly 
modified with respect to the non-irradiated one. Heating of the upper layers of 
the atmosphere leads to a nearly-isothermal external profile, yielding eventually 
a temperature inversion in the outermost layers^. As a consequence of this 
external heating, the radiative zone of the atmosphere extends to deeper layers, 
so that the radiative-convective boundary retreats towards the interior. This 
demonstrates the necessity to do consistent calculations between the irradiated 
atmosphere profile and the inner profile. Matching arbitrarily the inner profile 
to an atmospheric profile defined by the equilibrium temperature of the planet, 
as done in previous evolutionary calculations, yields severely flawed results. 



^Notc that the outermost temperature inversion observed in Jupiter, which afTects the UV part 
of the spectrum, is due to chromospheric cfTects, not included in the present calculations. These 
upper atmosphere effects, however, do not bear consequences on the evolution. 
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Figure 2. Spectrum of a planet with intrinsic temperature Tgfj = 500 
K, orbiting a G2 star at 0.3, 0.5 and 1.0 AU, respectively, from top to 
bottom, from Barman et al. (2001). The bottom curve displays the 
non-irradiated spectrum of the planet. These calculations correspond 
to a cloudless COND atmosphere model 



Figure 2 displays the emergent spectrum of a Tgff = 500 K, \ogg = 3.5 
planet orbiting at 0.3, 0.5 and 1.0 AU from its G2 parent star. The infrared 
part of the spectrum displays the intrinsic, thermal spectrum of the planet and 
exhibits the absorption features mentioned in §3. These features, however, are 
less pronounced than in the non-irradiated case, as seen from the comparison 
with the bottom line in the figure, a consequence of the aforementioned flat 
atmosphere temperature profile. This figure displays the calculations done with 
the so-called COND models, where grains are supposed to have sunk below the 
photosphere. This case maximizes the effect of irradiation, for the impending 
radiation penetrates deeper into the atmosphere than in the case of the so- 
called DUSTY models, where cloud layers form above the photosphere (see 
Barman et al. 2001). The short wavelength part of the spectrum corresponds 
to the reflected light from the parent star and clearly mirrors the absorption 
features of this latter. At short wavelengths, the albedo is dominated by H or 
H2 and He Rayleigh scattering. Because of the "grey" characteristics of dust 
opacity, the reflected spectrum in the case of the DUSTY models displays fewer 
spectral features (Barman et al. 2001). Detectability limits of the SIRTF and 
JWST missions are indicated in the figure for illustration. Note also the contrast 
between the planet and the star flux, which goes from ~ 10"'^ in the IR to larger 
values shortward of ~4 /im. 

Note that the present calculations assume that all the incident flux is de- 
posited at the substellar point, namely the planet point closest to the star, and 
that the energy is concentrated on the day-side of the planet. Such assumptions 
maximise the effect of irradiation. 
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5.2. Evolution 

First consistent evolution calculations of irradiated planets have been conducted 
by Baraffe et al. (2003) and more recently by Burrows et al. (2003). The 
calculations proceed as described earlier for the non-irradiated case. A grid of 
irradiated atmospheric structures is computed for various (log (7, T^s) conditions, 
for a given orbital distance a, and interpolation between these structures yield 
the unique boundary condition with the internal structure. As mentioned above, 
irradiation pushes the radiative-boundary limit further deep in the interior, so 
that, under strong irradiation conditions, the boundaty condition between the 
atmospheric and internal structure may lie in a radiative layer. For this reason, 
we have computed the Rosseland mean of the atmospheric opacities and used 
this value to calculate the radiative gradient for the interior structure. The 
incoming radiative flux Fine = l/2(-R*/a)F* modifies the energy budget of the 
planet and the emergent total flux Fout, which now reads (Baraffe et al. 2003): 

Font = Fine + crjj = A Fine + (1 " ^) -^inc + C T^S (3) 

where A denotes the Bond Albedo and the first term on the r.h.s. of (3) describes 

the reflected part of the spectrum. Wc stress that the albedo A is not a free 
parameter but is calculated consistently with the radiative transfer equation. 
This yields for the total luminosity: 

f dS 

Ltot = Lreflected + 47ri?^C7(Tg^^ -|- T^q) = Lreflected + '^'^Rl'^T^q - j ^"^^"^ (^) 

Note that only the last term on the right hand side of eq.(4) concerns the 

evolution of the intrinsic luminosity of the planet Lint- The second term defines 
its equilibrium temperature, i.e. the temperature the planet would reach in the 
absence of any internal source of energy, including from its contraction. The first 
term, which defines the Bond albedo, illustrates the luminosity fraction of the 
parent star reflected by the planet atmosphere. The dominant contribution to 
Lint, both in the irradiated and non-irradiated cases, arises from the contraction 
work contribution, ^^dm = - jj^^^dm. 

Figure 3 illustrates the evolution of the radius for different planets orbiting a 
G2 star, including the two detected transit planets HD209458B and OGLE-TR- 
56B. At a given orbital distance, namely a =0.05 AU, the effect of irradiation is 
smaller for the more massive, and thus hotter 10 Mj„p planet, compared with the 
0.69 Mjup one (HD209458B). This was expected from the smaller modification 
of the inner atmosphere profile due to a smaller stellar /planet flux ratio (see 
Fig. 5 of Baraffe et al. 2003). Irradiation, however, slows down substantially 
the contraction of the less massive planet, in particular at the early stages of 
evolution, when degeneracy effects in the interior are smaller. The flgure displays 
also the evolution of OGLE-TR-56B, for its recently revised mass, m ~ 1.45 it 
0.23 Mjup, for a radius 1.23 ± 0.15Rjup (Torres et al. 2003) and an orbital 
distance ~ 0.023 AU. The observed radii of the two transits are shown in the 
figure. The theoretical radius lies clearly below the observed one for HD209458B, 
but fits within the error bar for OGLE-TR-56B for the age of the system. We 
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Figure 3. Evolution of non- irradiated (dash) and irradiated (solid) 
gaseous planets orbiting a G2 star with Tgjj^ = 6000 K. The 0.69Mj>jp 
and lOMjup planets are located at 0.05 AU, while the LSMj^p is at 
0.023 AU 



derive m ~ 1.05 Mq, t ~ 4 Gyr to reproduce the observed extinction-corrected 
magnitude of the parent star OGLE-TR-56, in agreement with Sasselov (2003). 
These results show that there is presently no inconsistency between theory and 
observation for this transit (assuming OGLE-TR-56B is indeed a real transit). 

As discussed at length in Baraffe et al. (2003), however, no consistent evo- 
lution calculation can reproduce the observation of HD209458B within the error 
bars. Note that our calculations take into account the extension of the atmo- 
sphere, which we find to be of the order of ~ 0.05 Rj^p (at optical wavelengths) 
for the age of the object. Given the remaining uncertainties in the opacities 
of these planets, which should include processes such as non-equilibrium chem- 
istry, dynamics of grains, etc., we cannot exclude larger extensions. However, 
as mentioned above, the present calculations maximize the effect of irradia- 
tion, and it is very unlikely that different opacities modify the internal adiabat 
strongly enough to get a further 30% increase of the planet radius. Such a mod- 
ification requires a dramatic change of the energy deposited at the top of the 
internal adiabat, which determines the radius. But either molecular absorption 
or, at hotter temperatures, bound-free transitions and free-free scattering of the 
electrons in the interior prevent photons to penetrate deep enough. 

As discussed in §4.3 of Baraffe et al. (2003) and in Guillot & Showman 
(2002), one needs an extra source of energy to be deposited early enough in the 
planet, in the first ~ lO'^ years, and to remain constant after this stage, in order 
to reproduce the radius of HD209458B. Tidal dissipation or synchronization are 
unlikely, for these timescales are at most of the order of 10*^ yr. Note in passing 
that Bodenheimer et al. (2003) point out the possibility of a second planet 
orbiting HD209458B. This constant tidal heating could provide the required 
mysterious extra source of energy. 
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6. Brown dwarfs versus planets 

The deuterium-burning minimum mass, musMM — 0.12 Mq (Saumon et al. 
1996, Chabrier et al. 2000b) has been designated recently by the lAU as the 
official limit between BDs and planets. Although any definition can be accepted 
on a purely semantic ground, it is important to recall why the DBMM is sup- 
posed to play a particular role in star formation, distinguishing eventually two 
different formation mechanisms. This issue, as various star formation scenarios, 
is discussed at length in the review by Chabrier (2003) and will be only shortly 
summarized below. 

The particular role of the DBMM comes from the canonical star formation 
theory of Shu, Adams & Lizano (1987). In this theory, low-mass stars form 
from the collapse of initially hydrostatic, unstable dense cloud cores which have 
reached a p(r) oc r^^ density distribution of a singular isothermal spheroid. 
In this scenario, D-burning at pre-main-sequence ages is the key ingredient to 
trigger star formation. The onset of D-burning induces convective instability in 
the interior. Combined with the rapid rotation resulting from the accretion of 
angular momentum with mass, this convection is believed to generate a strong 
magnetic field through dynamo process. The field will drive a magnetocentrifu- 
gal wind that ultimately sweeps away the surrounding accreting material, and 
determines the mass of the nascent star. Within this picture, objects below the 
D-burning limit can not reverse the infall and thus can not form gravitationally 
bound objects. In this picture, the ambipolar diffusion timescale. Tad ^ 10 Myr 
for a cloud of average density ~ 10^ cm~^ and B a. few /itG, sets up the char- 
acteristic time scale for star formation. Such a scenario can now be reasonably 
excluded for several reasons: (1) SSOs are fully convective, with or without D- 
burning (Chabrier & Baraffe 2000, Chabrier et al. 2000a), (2) star formation is 
a rapid process, more rapid than the thermal crossing timescale or the aforemen- 
tioned ambipolar diffusion timescale, and comparable to the dynamical timescale 
Tdyn = (37r/32G/o)^/^ Ril-5xl0^ yr for typical star-forming molecular clouds, (3) 
in the ambipolar diffusion scenario, the cloud must survive long enough in near 
equilibrium between magnetic and gravitational pressure. This is not consis- 
tent with observations of rapid star formation and cloud dissipation, and with 
the observed turbulent nature of clouds. Indeed, equipartition between kinetic, 
gravitational and magnetic energy fails to reproduce the observed properties of 
molecular clouds, dominated by super-Alfvenic and supersonic motions, where 
kinetic energy dominates magnetic energy, with a decay timescale approximately 
equal to a dynamical timescale. Moreover, ambipolar diffusion models require 
large static magnetic field strengths (~30-100 //G) exceeding Zeeman estimates 
for low-mass dense cores (see e.g. Padoan & Nordlund 1999 and references 
therein). This indicates that, although magnetic field probably plays some role 
in star formation, it is unlikely to be a dominant process. On these grounds, 
D-burning seems to be excluded as a peculiar process in star formation. This 
is supported obscrvationally by the numerous free-floating objects which have 
been detected below the DBMM (Zapatero Osorio et al. 2000), although one 
must remain cautious about the age determination of the cluster and the related 
mass determination from the theoretical models (see Baraffe et al. 2002). In 
conclusion, the choice of the DBMM as the mass limit to distinguish BDs from 
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planets does not appear to be supported by physical considerations and must 
be considered as a purely semantic definition. 

A more physically-grounded minimum mass for star formation is the opacity- 
limit mass for fragmentation. This defines the theoretical minimum mass for 
cloud collapse between the initial, low-density phase, where compressional heat- 
ing is balanced by molecular line emission, with P oc p, and the subsequent, 
high-density (pc > lO^^gcm"^ for a cloud temperature T ^ 10 K) adiabatic phase 

where the gas becomes optically thick to infrared radiation, with P oc p^/^ (Lar- 
son 1969). This sets a minimum mass for fragmentation of the order of Jupiter 
mass (Boss 2001 and references therein). Revealing observationally the existence 
of this limit would be a tremendous breakthrough. 

The triggering process for star formation seems to be the small scale dissi- 
pation of large scale supersonic turbulence. This process leads to the formation 
of dense cores, some of which becoming dense enough for pressure forces to 
balance gravity (Padoan & Nordlund 2002). These cores eventually fragment 
further to form multiple systems (see e.g. Delgado-Donate et al. 2003, Goodwin 
et al. 2003). In terms of mass function, the mass distribution of the cores thus 
corresponds more or less to the IMF of systems, whereas further fragmentation 
leads to the IMF of individual objects. As shown by Chabrier (2003, figures 
3-5), the present observations of BDs either in the field or in young chistcrs are 
consistent (within possibly a factor of 2 to 3, depending on the real BD binary 
fraction) with the extension into the substellar regime of the same stellar single 
object or system IMF. In the opposite, the mass distribution of exoplanets differs 
drastically from such an IMF. A statistical analysis of their mass distribution, 
corrected for the uncertainty due to the inclination sin i of the orbital plane on 
the sky, has been established by Zucker Sz Mazeh (2001) and Jorissen, Mayor & 
Udry (2001). The resulting planetary mass distribution dN/dnip peaks around 
~ 1-2 Mjup - significantly below the DBMM - due to present detection limits of 
radial velocity surveys to detect smaller objects, and decreases rapidly with only 
a few objects above ^ 10 Mj„p. This distribution corresponds to a relatively flat 
MF, ^(logm) west, below ^ IQMjup. This different mass distribution clearly 
points to a different population of SSOs, namely planetary companions of stars, 
and thus to a formation mechanism different from the one which yields the stel- 
lar -|-BD IMF. The only identified difference between GPs like Jupiter or Saturn 
and BDs is the presence of a rock-l-ice core, of several Earth masses, at the very 
center of these planets. This ovcrsolar average abundance of heavy elements for 
our jovian planets suggests the presence of a disk in their formation environment 
and thus differs from the star-|-BD collapse scenario. Planets are thus unlikely 
to form as isolated, free floating objects, although some could have been ejected. 
On the other hand, the significant fraction of BDs in young clusters observed to 
harbor near and mid-IR excess consistent with dusty disks points to a forma- 
tion mechanism for BDs similar to the standard cloud collapse star formation 
process (Jayawardhana et al. 2003). The only indirect clue about the presence 
of a rocky core in our jovian planets, however, stems from the very accurate 
determinations of their gravitational moments from the Voyager, Pionneer and 
Galileo missions. Such data are obviously unavailable for the exoplanets (al- 
though transit observations may provide some information about the presence 
of a core) and for the moment the only possible way to distinguish BDs from 
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planets is a statistical argument, with two radically different mass distributions. 
Therefore, the observed so-called free floating "planets" or "sub brown dwarfs" 
represent most likely the low-mass tail distribution of usual brown dwarfs. In 
any case, the search for such very-low mass objects is of tremenduous interest 
for it will allow us to identify the minimum mass for star formation. 

7. Conclusion and perspective 

This review outUncs our present understanding of the physics of substcllar ob- 
jects, brown dwarfs and irradiated or non-irradiated gaseous planets. Various 
successes of the theory, which had correctly predicted different mechanical and 
thermal properties such as the EOS up to tlie mcgabar range, the spectral energy 
distribution of G1-229B like objects, dominated by methane absorption, the pho- 
tometric evolution of dusty and dustless atmosphere SSOs, and the mass-radius 
and mass-magnitude relationships in the stellar range, give us confidence in the 
reliability of its basic foundations, although there is certainly room for a lot 
of improvements. As argued in this review, a consistent treatment between the 
atmospheric profile and the interior profile is mandatory to yield reliable results, 
both in the irradiated and non-irradiated case. This theory allows us to derive 
the fundamental properties of an object (mass, temperature, age, radius) from 
its observable spectroscopic or photometric signatures. Conversely, the theory 
provides reliable basis to guide observational projects. Our understanding of 
the properties of strongly irradiated planets, the so-called hot jupiters, is still 
questioned. Whereas the theory adequately reproduces the radius of OGLE-TR- 
56B, no calculation can reproduce the one of HD209458B. Whether this latter is 
a common representative of hot-jupiter objects or whether its large radius is due 
to tidal interaction with an undetected companion needs to be assessed. Ongo- 
ing searches for transit planets will help to answer this question and confirm or 
not our present description of strongly irradiated planets. 
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